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Abstract. We report relationships between the effects of noise and ap- 
plied constant currents on the behavior of a system of excitable elements. 
The analytical approach based on the nonlinear Fokker-Planck equation 
of a mean-field model allows us to study the effects of noise without ap- 
proximations only by dealing with deterministic nonlinear dynamics . We 
find the similarity, with respect to the occurrence of oscillations involv- 
ing subcritical Hopf bifurcations, between the systems of an excitable 
element with applied constant currents and mean-field coupled excitable 
elements with noise. 
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1 Introduction 

The effects of noise on neural systems are one of the great concern for researchers 
in neurosciences. Among theoretical and computational studies is the influence 
of noise on excitable systems as well as sustained oscillatory ones. Simulations 
of coupled Hodgkin-Huxley neurons subjected to noise are one of typical exam- 
ples to demonstrate its behavior (T]. Numerical solutions of the Fokker-Planck 
equation of coupled active rotators with excitatory and inhibitory elements have 
been obtained and analysed [2]. Another approach of Gaussian approximations 
has also been used to reveal bifurcation structures with changes in noise inten- 
sity for investigating a system of coupled active rotators and Fitz-Hugh Nagumo 
neurons |3I4) . 

In this paper, we propose a model of coupled excitable elements under the 
influence of external Langevin noise on the basis of mean- field concept [51617] . 
Satisfying self-averaging property, the model enables us to describe the time 
evolution of order parameters of the system without approximation. Investigating 
this dynamic system instead of a set of Langevin equations, we show similarity 
between noise and applied currents effects on coupled and uncoupled excitable 
elements from the viewpoint of bifurcation structures. 
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2 Excitable Element 
2.1 Model 

First, we consider the effects of applied constant currents on a system of an 
excitable element. Let us suppose that dynamics of a 2-dimensional system of 
z ( x \ z^) satisfies the following set of differential equations: 

a (x) z (x) +J (x) F (x)( b (x,x) z (x) + b (x,y) z (y)j + j ^ 

a (v) z (y) + j(y)F^(b iy ' x) z ( - x) + b (y - y) z (y) ), (2) 

where b^ ,l/ \ (/i = x,y) are constants, / is the applied current and 
F^{-) are coupling functions. To make the system excitable, we specify F^ x \-) 
and F^ v \-) as nonlinear and linear functions as 

F^(z) = zexp(~Y (3) 

F^{z)=z. (4) 

The model parameter values are = 2.5, = 0.0030, = 1.5, b^ = 

0.50, b^^ = 4.0, b( y ^ = 1.0, J( x ) = 5.0, = -0.0040. In this situation, z<» 
and z( y > correspond to fast and slow variables and the system exhibits excitable 
properties for the applied constant current I = —3.0. 

Figure [1] shows the behavior of the system. When the system is perturbed, 
z (x) ig e ^ s exc it ec i' to generate a single pulse, while z^ exhibits periodic pulses 
with a certain level of applied constant currents / (Fig.Q](a)). We can also easily 
understand this dynamics in phase plane (Fig.[T](b)). Taking an initial condition 
right side from the equilibrium point, the trajectory begins its journey to return 
to the fixed point. 



dz<- x > 

dt 
dz^~> 

dt 



2.2 Bifurcation Diagrams against Applied Constant Currents 

To reveal the bifurcation structures with changes in applied constant currents /, 
we conduct a linear stability analysis (Fig. [5]). Let us start with the case / = —3.0, 
where the fixed point of the system is the stable node. With increasing the ap- 
plied constant current /, this point moves toward origin and its stability changes 
to the stable spiral, to the unstable spiral, and to the unstable node, yielding 
stable limit cycles accompanying the subcritical Hopf bifurcation. Note that the 
fixed points have rotational symmetry as (zq X \zq V \i) = (— z^ x \ — — /). In 
this way, appropriate levels of applied constant currents \I\ < 2.4038 give rise to 
oscillatory states, whereas these oscillations disappear in higher or lower ampli- 
tude |/| > 2.4042. In amplitude 2.4038 < \I\ < 2.4042, the limit cycle and fixed 
point attractors coexist. 
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Fig. 1. Dynamics of the system of an excitable element, (a) The time evolution 
of z( x >. gets excited once by the perturbation at t = 300.0. After the changes 
in amplitude of the applied constant current (f > 1000.0), begins to oscillate, 
(b) Nullclines in phase plane. At I — —3.0, the trajectory converges to the stable 
node after sufficiently large times. 



3 Coupled Excitable Elements under the Influence of 
Noise 

3.1 Model 

To compare the effects of noise and applied constant currents on excitable ele- 
ments, we extend the system described by Eqs. ((T|) - ^ to a coupled one. We 
suppose a system of N excitable elements and assume that they interact with 
each other via coupling functions FM^) with denoting coupling constants. 
In addition, each element of a system are subjected to the independent Langevin 
noise ^[^(t). Each of those systems reads 

= -aWzP + ±- V J^F^Xb^zf + b^zf) + I + 4 x \t), (5) 
at 1\ £ — ' J J 

j (v) i N 

= - a (y)z^ +±yj(y)F(y\b^ x h^ ) + b^yh ( j y) ) + 4 y) (t), (6) 

at N J J 

i=i 

where i = 1,---,N. We postulate that the Langevin noise i]^(t) are white 
Gaussian ones, {^(t)) = 0, {r^ {t)r}f ] (*')) = 2D^§ ij 5 llu 8{t - if). The cou- 
pling functions F^(-) and the model parameters are the same as described 
above. The applied constant current is set to / = —3.0. In the absence of noise, 
Eqs. © and © with N = 1 recover Eqs. JU and ©. 

3.2 Nonlinear Fokker-Planck Equation Approach 

In the thermodynamic limit N — >• oo, we can take advantage of reducing Eqs. ([3]) 
and (0 to a single body equation as seen below. The mean-field coupling terms 
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Fig. 2. Bifurcation diagrams of the system of an excitable element, (a) The 
whole diagram. H u and dc denote the subcritical Hopf bifurcation and double 
cycle, respectively, (b) Magnified view of the occurrence of the limit cycle. The 
limit cycle appears accompanying the subcritical Hopf bifurcation. 



satisfy the self-average property, which is written by the empirical probability 
density P(t, z<», z^), 



(J*))= / dz^dz^F^\b^' x h^ + b M z^)P(t,z^ x \z^). 



(7) 



Then the system of Eqs. ((3J) - ([BJ is indeed reduced to the 1-body dynamics 

Av) as 



dt 
dz^ 
dt 



= -a^zM +J (x) {F (x) ) +1 + C (B) (<), 
= -abhto +J {v) {F {v) ) + ( {y) (t), 



with white Gaussian noise C (/l) (*)» (C (/i) (*)) = °> (C M (*)C (l,) (*')> = 2D W V^- 
£'). Thus, one obtains the nonlinear Fokker-Planck equation for the empirical 
probability density corresponding to the above Langevin equations [5], 

d 



d 



d 



-a {x) z {x) + j( x )(F (x) ) + / - D {x) 
-a {v) z {v) + J {v) (F (y) ) ~ D {v) 



dziv) 



P. 



P 



(8) 



A Gaussian probability density is a special solution of the nonlinear Fokker- 
Planck equation © . Furthermore, since the H theorem [5] ensures that the prob- 
ability density satisfying Eq. © converges to the Gaussian-form for sufficiently 
large times, we are concerned with the Gaussian probability density as 
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P G (t,Z 



(z) r (j/h 



1 



: exp 



27r v /detC G (t) 



where (-) G denotes expectation over P G . Then, the coupling terms Eqs. d7J) are 
described only up to the second moments. We derive a set of closed ordinary 
differential equations as 



d{z^) G 

dt 
d(z(»>) G 

dt 

dt 

d(u^ 2 ) c 
dt 



-a^(z^) G + J^(F^) G + I, 
>)(zW) G + jW(F(")) G , 
-2a (a;) (u (:l;)2 } G + 2L> (:r) , 



d<u y ))a =-(aW+qW)<uWtiW) a> 



where one has from Eqs. © - ((U) 



V /G (<7 2 + l) 3 / 2 



exp 



2(a 2 + 1) 



(F 



(y)\ — miv) 



(9) 
(10) 

(11) 
(12) 

(13) 

(14) 
(15) 



and = b^ x \z^) G + b^{z^) G , a 2 = b (x ' x ^{u^ 2 ) G + b^ x ' yS > 2 (u^ 2 ) G . 

Note that (u^u^) G -> and (m (ai)2 ) g -> D^/a^ (t -> oo), implying that 
the external Langevin noise contributes to the dynamics through the variance. 
For simplicity, we assume the Langevin noise intensity £)W = in what follows. 

Investigating the dynamic system of the order parameter equations written 
by Eqs. (!§))- (fT?|) together with Eqs. (fl"4"f and (1T51) . we can compare the effects 
of noise and applied constant currents on the systems of globally coupled and 
uncoupled excitable elements. 

The dynamics of order parameters of the system is shown in Fig. [3] While 
the mean value (z^ x ') G converges to the stable fixed point for lower intensity of 
the Langevin noise D^ x \ it oscillates for relatively higher noise intensity 
(Fig. [3] (a)). This oscillatory state suggests that individual excitable elements 
almost simultaneously get excited under the influence of noise, i. e. noise-induced 
synchronization in coupled excitable elements. The portrait of phase space (Fig. 
[3] (b)) is similar to the case of uncoupled deterministic system (Fig. Q] (b)). 
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Fig. 3. Dynamics of the system of coupled excitable elements under the influence 
of noise, (a) The time evolution of (z^ x ') G . Whereas z^ x ' does not oscillate for 
lower noise intensity D( x > (t < 500.0), it begins to oscillate for an appropriate 
noise level (t > 500.0). (b) Nullclines in phase plane with the noise intensity 
£)i x ) = 0.4. In both figures, the applied constant current is set to / = —3.0. 



3.3 Nonequilibrium Phase Transitions and Synchronization Induced 
by Langevin Noise 

Using the nonlinear Fokker-Planck equation approach with the mean- field model, 
we can identify the occurrence of genuine bifurcations of order parameters with 
changes in noise intensity. The result of numerical survey is shown in Fig. 01 

In the deterministic limit D( x > — »• 0, the fixed point of the system is a stable 
node. With increasing D^ x \ this stability changes to a stable spiral and further 
to an unstable spiral. A limit cycle attractor occurs accompanying the subcritical 
Hopf bifurcation (Fig. [4] (b)), similar to the deterministic uncoupled case with 
changes in applied constant currents (Fig. [5] (b)). In contrast, the disappearance 
of the oscillatory state is the result of the supercritical Hopf bifurcation (Fig. 
H] (c)). So it is clearly seen that a certain level of noise induces synchronized 
oscillations in the system of mean-field coupled excitable elements under the 
influence of noise. 



4 Summary 

We have studied relationships between the roles of noise and applied constant 
currents in a system of excitable elements. The analytical approach using the 
nonlinear Fokker-Planck equation associated with the mean-field model allows us 
to obtain genuine bifurcation diagrams against noise intensity. In a system of an 
excitable element without noise, applied constant currents of a certain amplitude 
have given rise to oscillatory states, which originate from the subcritical Hopf 
bifurcation. The similar structure has been found in coupled excitable elements 
under the influence of noise. 
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Fig. 4. Bifurcation diagrams of the system of coupled excitable elements un- 
der the influence of noise, (a) The whole diagram. H s and H u , and dc denote 
the supercritical and subcritical Hopf bifurcation, and double cycle, respectively, 
(b) Magnified view of the occurrence of the limit cycle. The limit cycle appears 
accompanying the subcritical Hopf bifurcation, (c) Magnified view of the disap- 
pearance of the limit cycle. The limit cycle vanishes with the supercritical Hopf 
bifurcation. 



Further analysis of the effects of the multiplicative noise on a system of 
mean-field coupled excitable elements will be reported elsewhere. 
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